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Definition 1. Algorithms which change the boundary of the solution region in order to find 
the optimal solution of an integer programme are called cut algorithms. The branch-and-bound 
algorithm does this by splitting the solution region into two and then discard the one which does 
not contain the optimal solution. The Gomory algorithm, on the other hand, reduces the feasible 
region with the help of a new constraint without the region being splitted. 
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Definition 2. We call branching a process by which a programme whose solution contains a 
non-integral 7 < x; < k is made into two separate programmes having the additional constraint 
x; <j inone, and x; > k in the other, the objective together with all the constraints of the original 
problem of which remain the same. Here j and k are positive integers and j < k. 
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Definition 3. In the branch-and-bound algorithm, if the objective is maximisation, the value of 
the objective obtained when the first integral approximation occurs is said to be the lower bound 
for the problem, and if the objective is minimisation it is said to be the upper bound of the same. 
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Algorithm 1 Branch-and-bound algorithm for integer programming. 


find first approximation 
while approximations not all integers do 
choose x; from all non-integral variables such that min (|x; — |a;||, |e; — [a:]|) 


is maximised 


branch 
choose the branch whose value of the objective is maximum 


endwhile 
solution<last approximation 


Example 1. (Problem 6.9; Bronson, 1982) 


maximise: z= 2, +2724 73 
subject to: 2%, +3224 323 < 11 
with: all variables non-negative and integral 


Solve by branch-and-bound algorithm. 


Solution. Draw a simplex table of Programme 1. 


Ty v2 v3 v4 
1 2 0 
r4 0 2 3 3 1 11 
-1 -2 -1 0 0 
Replace x4 for x2 as the basic variable. 
Ly x2 x3 LA 
Ps a x 
iv a 3 z 
l= 4 = 3.6, ec} = x3 = x4 = 0, 2* = Since 3 < x3 < 4, branch into two programmes, 


namely Programme 1 where x2 < 3, and Programme 2 where x2 > 4. Consider first Programme 2. 


maximise: z=2, +2724 73 
subject to: 2%, +3294 323 < 11 
r2 <3 
with: all variables non-negative and integral 
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Use the simplex method in a tabulated form. 


Ty r2 v3 v4 v5 
1 2 1 0 0 
L4 0 2 3 3 1 0 11 
X5 0 0 1 0 0 1 3 
-1 -2 -1 0 0 0 
Replace the basic variable x5 with x2. 
Ly x2 x3 L4 X5 
L4 2 0 3 1 -3 2 
XQ 0 1 0 0 1 3 
-1 0 -1 0 2 6 
Replace the basic variable x4 with 7x1. 
Baal v2 v3 TA x5 
1 1 o 8 2? 8 1 
XQ 0 1 0 0 1 3 
0 0 z 4 s 7 
ej =1, 45 = 3, 73 = 2, = 25 = 0, z* =7 Then consider Programme 3. 
maximise: z= 2,1 +2724 73 
subject to: 22, +3224 323 < 11 
w>4 
with: all variables non-negative and integral 
Draw a table for the two-phase method. 
ZY r2 v3 v4 v5 a) 
1 2 1 0 0 —M 
L4 0 2 3 3 1 0 0 11 
XE —M 0 1 0 0 -1 1 4 
-1 -2 -1 0 0 0 0 
0 -1 0 0 1 -1 -15 
Change x4 for x2 in the basic variables. 
Ly x2 x3 4 £5 x6 
r2 2 1 1 3 0 0 + 
x6 — 2 0 -1 — 3 -1 1 5 
5 0 1 3 0 0 2 
3 0 1 $ i a — 


The coefficient parts of the row corresponding to the basic variable xg and the last row cancel 
each other. The optimal result is 73 = oan Ly = 23 = 24 = LE = 2G =O and z* = 2. 
(3) LQ 2 4, a= 2, (0, 3) 


(2) a2 <3, 2* =7, (1,3) 


* 


Therefore the solution is ef = 1, x3 = 3, x3 = 7] = x5 = 0, and z* = 7. 


Algorithm 2 Gomory algorithm for integer programming. 


while solution not wholly all integers do 
choose one non-integral optimal approximation 
write a relation from the row where that variable is basic 
rewrite the relation to make all fractional coefficients some integer plus a proper fraction 
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move all the fractions to LHS, and all the non-fractions to RHS 

write a new constraint as LHS > 0 

find the solution for the original problem together with the new constraint 
endwhile 


Example 2. (Problem 7.1; Bronson, 1982) 
maximise: z= 27,4 22 
subject to: 22, + 5% <17 

321 + 2% < 10 


with: 21,22 non-negative and integral 


Use cut algorithm. 


Solve 


Solution. Find the first approximation of Programme 1 normally using the simplex method. 


Ty r2 v3 T4 
2 1 0 0 
x3 0 2 5 1 0 17 
vw, 0 3 2 0 1 10 
-2 -1 0 0 0 


Since 2 < zZ, we know that 3 is the pivot element, and therefore we replace the basic variable 
v4 with 71. 


1 v2 v3 v4 
o a a z 
a a re 8 
0 3 O 3 3 
We have x7 = 2, x3 = a, v5 = 24 =O and z* = 2. Since both x} and x3 are non-integers, 
arbitrarily choose the former to generate a new constraint. Then our Programme 2 becomes, 
2 1 10 1 
+ 3% + 384 = 3 = 34 3 
2 1 1 
3% 4 34 a 3-2) 
2 1 1 
gtr atte 
1 1 
gta hats 
2% +44 >1 


and our new programme becomes 
maximise: z= 27, + 272+ 0273+ 0x4 
11 2 3l 


subject to: = 
subject to 3 2 34 3 
25 _ 10 
Ly 32 374 = 3 
2%2 +44 >1 


with: all variables non-negative and integral 


ry x2 w3 v4 U5 v6 
2 1 0 0 0 —-M 
Ly 0 1 2 0 3 0 0 2% 
r3 0 0 y 1 —2 0 0 4 
r6 —-M 0 2 0 1 -1 1 1 
-2 -1 0 0 0 0 0 
0 -2 0 -1 1 -1 -1 
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Now 22 replaces zg in the basic variables and becomes the pivot element. 


1 v2 v3 v4 5 U6 
v1 1 0 0 $ —3 3 
a3 0 1 is ui i 1 
ra ee ee ! 
2 0 0 4¢ -- 4 ! 
0 0 0 0 0 0 0 
This becomes, 
Ty v2 v3 LA v5 
X1 1 0 0 0 3 3 
x3 0 0 1 —3 ¥ ¥ 
tty 0 1 0 2 Ak u 
a oe 2 
Then our first approximation of Programme 2 is «7 = 3, £3 $ xe iu, xv, = «5 = 0, and 
= B. Arbitrarily choose x3 to generate the new constraint. 
ie 1 Le at 
v2 gv gs = 5) 
1 1 1 
On te 
1 1 1 Sa 
Ye es = 
94 95 a= 
t4—-2%5>1 
Then our Programme 3 becomes, 
maximise: z= 2%, +2%2+02%3 + 0x4 + Ov5 
1 
subject to: 2, + 3% =3 
rr Thy AG 
mee a toa, a 
1 1 
Pa Gta GS = 5 
t4—-%5>1 
with: all variables non-negative and integral 
We draw our table for this programme. 
1 v2 v3 v4 U5 v% T7 
2 1 0 0 0 0 —M 
v1 0 1 0 0 0 5 0 0 3 
x2 0 0 1 0 4 —4 0 0 $ 
x3 0 1 0 0 -3 i 0 0 u 
x7 —M 0 0 0 1 -1 -1 1 1 
-2 -1 0 0 0 0 0 0 
0 0 0 -1 1 1 -1 -1 


Then 2x4 replaces the basic 77 to become the pivot element. 
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Ly r2 13 ra w5 v6 
ry 1 0 0 0 3 0 3 
rp 0 1 0 0 0 3 0 
v3 1 0 0 0 —2 —3 11 
r4 0 0 0 1 -l - 
-2 -1 0 0 0 0 0 
Next, x, remains basic and becomes a pivot element. 
T1 £2 v3 v4 v5 v 
L1 1 0 0 0 ; 0 3 
re 0 1 0 0 5 0 
x3 0 0 0 0 -1 -3 8 
L4 0 0 0 1 -1 -1 1 
0 -1 0 0 3 0 6 
This becomes 
L1 XQ v3 v4 Ls x6 
L1 1 0 0 0 3 0 3 
re 0 1 0 0 5 0 
x3 0 0 0 0 -1 —3 8 
L4 0 0 0 1 -1 -1 1 
0 0 0 0 3 $ 6 


The optimum point for Programme 3 is then, ej = 3, 73 = 8, 7] = 1, 7 = #5 = 7% =0 
and z* = 6. Therefore the solution to the original problem Programme 1 is «jf = 3, x5 =0 at the 
objective value z* = 6. 


# 
Definition 4. A transportation problem involves m sources each of which supplies a;, i = 1,...,m, 
units of a certain product, and n destinations each of which requires b;, i = 1,...,n, units of the 


same. The problem may be stated as following. 


m nm 
maximise: 25 ) Cig Lig 


i=1 j=1 


n 
subject to: y Lig=a, t=1,...,m 
j=l 


m 
yoy Sx Yee reef) 
i=1 


with: all 2,; non-negative and integral 


The total supply and the total demand are assumed to be equal. Were this not so, a fictitious 
destination or a fictitious source is added. 
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Definition 5. The north-west corner rule finds an initial basic solution for the transportation 
algorithm of the integer programming. It begins with the (1,1) cell in the m x n table, and allocates 
as many units as possible to x1, violating neither the constraints of supply, that is the summation 
along each row, nor those of demand, that is the summation along each column. Then carry on 
moving for each step either right or downwards, until we reach the lower-right corner, &nn. 
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Definition 6. A loop, which is a sequence of cells in the table used for finding the solution in the 
transportation problem, has the following properties. 
a. each pair of consecutive cells is on either the same row or the same column 
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b. no three, or in fact any odd-numbered, consecutive cells lie in the same row or column 
c. the first and the last cells are on the same row or column 
d. the path along the loop is self-avoiding, that is no cells appear more than once in the sequence 
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Algorithm 3 Transportation algorithm. 


while optimal solution not attained do 
find an initial, basic feasible solution using, for instance, the North-west corner rule 
let either uj = 0 or v; =0 depending on whether the i**-row or the j*®-column 
has the maximum number of basic solutions 
find all u; and v;,2=1,...,m andj = 1,...,n from u; + vj = cj; for basic variables, 
and from ¢;; — u; — v; for non-basic variables 
improve the solution 
endwhile 


Note 1. Ina transportation problem, optimal solution is achieved when c;; — ui; — u; > 0 for all 
transportation costs per unit c,; of all non-basic variables. 
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